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We perform a measurement of the CP asymmetry in b — ► S7 decays using a sample of 383 x 10 
BB events collected by the BABAR detector at the PEP-II asymmetric B factory. We reconstruct 
sixteen flavor-specific B decay modes containing a high-energy photon and a hadronic system X s 
containing an s quark. We measure the CP asymmetry to be —0.011 ± 0.030(stat) ± 0.014(syst) for 
a hadronic system mass between 0.6 and 2.8 GeV/c 2 . 
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The decay b — + S7 is a flavor-changing neutral cur- 
rent process described by a radiative penguin diagram in 
the Standard Model (SM). It is sensitive to new physics 
which can appear in branching fraction or CP asymmetry 
measurements. Measurements of the branching fraction 
[3, [|J are in good agreement with the SM Q predictions. 

A CP asymmetry between b — > and b — > 37 de- 
cays is predicted by the SM to be < 1% [J] but could be 
enhanced up to 15% [!, H, 0] in models of physics be- 
yond the SM. Existing measurements are consistent with 
zero CP asymmetry with a precision of 5% The in- 

creased precision obtained in this work allows us to better 
discriminate between various theoretical models 

We use a sample of 383 x 10 6 BB pairs collected at 
the T{AS) resonance by the BABAR detector [11| at the 
PEP-II e + e~ B factory. In addition, we use 36.3 fb _1 
collected 40 MeV below the T(AS) resonance to study 
backgrounds from non-i? decays. 

We reconstruct 16 exclusive b — ► S7 final states: 

B~ -> K s TT--/,K-n Q -/,K-7r + 7r--f,K s 7r-TT -f, 
iW7r°7r°7, K®ir + TT~Tr~j, iW7r + 7r~7r°7, 
i^7r-7r°7r°7, K~r]j, K + K~K-j, 

B° -> K-Tr + -f,K-Tr + TT "/,K-n + TT-TT + "f,K-Tr + TT 7r "f, 
K~n+r]^, K + K-K-ix + ll 

and measure the yield asymmetry with respect to their 
charge conjugate decays b — > S7 . These modes are se- 
lected because the particles in the final state identify the 
flavor of the B meson and they can be reconstructed with 
high statistical significance. 

The high-energy photon from the B decay is recon- 
structed from an isolated energy cluster in the calorime- 
ter, with a shape consistent with the electromagnetic 
shower produced by a single photon, and an energy 
E* > 1.6 GeV in the T(4S) center-of-mass (CM) frame. 

The hadronic system X s , formed from the kaons and 
pions, is required to have an invariant mass Mx s between 
0.6 and 2.8 GeV/c 2 , corresponding to a photon energy 
threshold E~ > 1.9 GeV in the B meson rest frame. 

Charged kaons are identified by combining informa- 
tion from the Cherenkov detector and the energy-loss 
measurements from the tracking system. The remaining 
tracks are assumed to be charged pions. The K® can- 
didates are reconstructed by combining two oppositely 
charged pions with an invariant mass within 9 MeV/c 2 of 
the nominal K® mass fl2| and a minimum flight distance 
of 2 mm from the primary event vertex. Both charged 
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and neutral kaons are required to have laboratory mo- 
menta > 0.8GeV/c. 

Neutral pions and rj candidates are reconstructed from 
pairs of photons with energies above 50 MeV in the labo- 
ratory frame and a lateral moment [l3| less than 0.8. The 
lateral moment measures the spread of a shower in the 
calorimeter and provides good separation between elec- 
tromagnetic and hadronic showers. The invariant mass 
of the pair of photons is required to be between 115 and 
150 MeV/c 2 for tt° candidates and between 470 and 620 
MeV/c 2 for 77 candidates 

Monte Carlo (MC) samples based on EvtGen [ll| 
and GEANT4 hj are used to simulate the signal and 
background processes and the detector response. The 
6 — * S7 signal sample is generated with a photon spec- 
trum derived from Ref. [j| assuming to& = 4.65GeV/c 2 . 
The fragmentation of the X s system is modeled using 
JETSETplI corrected to fit the BaBAR data as described 
later. 

The background to the B reconstruction is dominated 
by continuum processes (e + e~ — > qq, with q — u,d, s, c) 
that produce a high-energy photon either by initial-state 
radiation or from the decay of ir° and 77 mesons. Contin- 
uum events tend to be less isotropic than i?-decay events 
since they result from hadronic fragmentation of high- 
momentum quarks back-to-back in the CM frame. High- 
energy photons in these events tend to be collinear with 
the thrust axis formed from the rest of the event (ROE), 
defined as those particles not used in reconstructing the 
signal B candidate. We reject such backgrounds by re- 
quiring that the cosine of the angle between the photon 
and the thrust axis of the ROE (in the CM frame) be 
less than 0.85. We further reject the continuum events 
by requiring the ratio of the second (L2) and zeroth (Lo) 
Legendre moments for the ROE particles with respect to 
the B flight direction to be smaller than 0.46. 

Continuum events with high-energy photons from ir° 
and 77 decays are major backgrounds. To veto these 
events, we associate each high-energy photon candidate 7 
with another photon candidate 7' in the event. For mul- 
tiple 7' candidate in an event, we choose the 77' pairs 
whose invariant mass, determined from adding the four 
vectors, is closest to the nominal ir° mass (or 77 mass in 
case of ry veto). Events are rejected if the photon pairs 
are consistent with ir° or 77 decays based on the output 
of a boosted decision tree (BDT) [ljj constructed from 
the energy of the less energetic photon 7' and m 77 /. 

We reject the remaining continuum events by con- 
structing an additional BDT that combines information 
from a number of variables related to the event shape, 
the kinematic properties of the B meson, and the flavor- 
tagging [181 ] properties of the other B meson in the event. 
Examples of these variables are the Fox- Wolfram mo- 
ments [lH, and the cosine of the B flight direction com- 
puted in the CM frame with respect to the beam axis. 
Optimization of the selection criteria of the ir° veto, 77 
veto, and event selection BDTs is performed using an 
iterative method which maximizes the statistical signal 
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significance. After the final event selection, we reject 97% 
of the continuum background while retaining 55% of the 
signal events. 

Fully reconstructed b — > sj decays are character- 
ized by two kinematic variables: the beam-energy sub- 
stituted mass toes = y s/4 — p% 2 , and the energy dif- 
ference between the B candidate and the beam energy 
AE = E* B — \fs/2, where Eg and p* B are the energy and 
momentum of the B candidate in the e + e~ CM frame, 
and y/s is the total CM frame energy. Signal events are 
expected to have a AE distribution centered near zero 
and a toes distribution centered at the mass of the B 
meson. For events with multiple B candidates, we select 
the one with the smallest \AE\. 

We perform a one-dimensional fit of Toes to the data 
in the entire Mx a region ([0.6, 2.8] GeV/c 2 ) as well as in 
five different regions of M Xs ([0.6, 1.1], [1.1, 1.5], [1.5, 2.0] 
and [2.0, 2.8] GeV/c 2 ) to study whether the asymmetry 
has significant mass dependence. Only candidates in the 
range \AE\ < 0.10 GeV and 5.22 < to ES < 5.29GeV/c 2 
are considered. Probability density functions (PDFs) are 
constructed for both signal and background in the five 
Mx s regions. We use the charge of the reconstructed 
final state (B~/B + ) or the charge of the kaon (B°/B°) to 
define two flavor categories, and perform a simultaneous 
fit for the flavor asymmetry in each Mx 3 region. 

The signal events are described by a function 

/(toes) = exp[-(rn E s-Mo) 2 /(2cr| :fl , + a L ,_ R (TOES-Mo) 2 )] 
where the parameters are determined by an unbinned fit 
to the signal MC. In the above function, /io is the peak 
position of the distribution, ol,r are the widths on the 
left and right of the peak, and cxl.r parameterize the tail 
on the left and right of the peak, respectively. 

The background surviving the final selection can be 
attributed to one of three sources: continuum events, 
BB events other than b — > sj decays (referred to as 
generic BB), and "cross-feed events", defined as events 
containing a b — > S7 decay, but in which the true decay 
was not correctly reconstructed. The shape of the cross- 
feed and BB background is described by a binned PDF, 
determined from MC with 1 MeV/ c 2 binning. 

The continuum background is described by an ARGUS 
function [2(| determined from a fit to the off-resonance 
data. In this fit, the toes distribution is shifted to have 
the same end-point as that of the on-resonance data. 

In the maximum-likelihood fit, all parameters are fixed 
with the exception of the normalizations of the various 
components as well as fio, which is determined from fit- 
ting the data, since the peak position is not well modeled 
in the MC simulation. The signal, BB and cross-feed 
shapes are constrained by the MC, while the continuum 
background shape is fixed to that of off-resonance data. 
The shapes of the distributions are assumed to be the 
same for B and B candidates, with the exception of the 
BB and cross-feed background, which are allowed to vary 
between b and b in order to eliminate the possibility of a 
false CP asymmetry. In Figure [1] we present the final fits 
to the Toes distributions for b — > sj and b — > events 
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FIG. 1: Fits to the tties distribution in data for b — > S7 events 
in M x „ region (a) [0.6, 1.1], (b) [1.1, 1.5], (c) [1.5, 2.0], (d) 
[2.0, 2.8], and b -> sy events in M Xs region (e) [0.6, 1.1], 
(f) [1.1, 1.5], (g) [1.5, 2.0], (h) [2.0, 2.8], The dashed line 
shows the shape of the continuum, dotted-dashed line shows 
the fitted signal shape and the dotted line shows the BB and 
cross-feed shape. 



for the four Mx s sub-regions. As expected, the signal 
to background ratio decreases from lower to higher Mx s 
regions. In Figure [H we present the final fits to the to-es 
distribution for b — > S7 and b — > S7 events for the entire 
Mx B region. 

The direct CP asymmetry is calculated as 



A C p 



1 

CD) 



Nb + Nj- 



- AD) - A 



del 



(i) 



and 



where Nb and Nt are the yields of the b — > 
b — > S7 signals respectively. Adet, described in details 
below, is the flavor bias caused by the detector responses 
to positively and negatively charged particles. Table U 
presents the fitted values for (N b - N^)/(N b + N%). 
AD = ([lu — to) is the difference in the wrong-flavor 
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Mxs 

( Gev/ c ) 


iV 6 -JV- 
AT . 4- AT— 


Adet 


BB and cross-feed 
model syst 


Continuum 
model syst 




0.6-1.1 


O.Olo ± 0.029 


O.OOo ± 0.014 


0.002 


n nn 1 
0.004 


0.010 ± 0.029 ± O.Olo 


1.1-1.5 


-0.003 ± 0.049 


-0.003 ±0.015 


0.003 


0.004 


0.000 ± 0.049 ± 0.016 


1.5-2.0 


-0.064 ±0.077 


-0.017 ±0.010 


0.010 


0.002 


-0.047 ±0.077 ±0.014 


2.0-2.8 


-0.097 ±0.180 


-0.002 ±0.005 


0.070 


0.168 


-0.077 ±0.180 ±0.182 


0.6-2.8 


-0.018 ±0.030 


-0.007 ±0.005 


0.012 


0.006 


-0.011 ±0.030 ±0.014 



TABLE I: For each Mx B bin, we present the fitted CP asymmetry: (Nb — N-j-)/(Nb + Nr), the flavor-bias of the detector: 
Adet, the systematic error arising from the BB and cross-feed modeling and the systematic error arising from the continuum 
background modeling. The last column shows the final results for the CP asymmetries. 




FIG. 2: Fits to the tties distribution in data for (a) b — > s'y 
events and (b) b — > 57 events in the entire Mx s region. The 
dashed line shows the shape of the continuum, the dotted- 
dashed line shows the fitted signal shape, and the dotted line 
shows the BB and cross-feed shape. 



fraction between b and b decays, and (D) = 1 — (p + 
lo) is the dilution factor from the average wrong-flavor 
fraction. The small wrong-flavor fraction Co (lo), defined 
to be the fraction of b (b) reconstructed as the opposite 
flavor, is due to charged pions misidentified as charged 
kaons. Using the particle misidentification rate measured 
in control samples in data we calculate AD = (5 ± 4) x 
IfJ- 5 and 1 — (D) = (5.4 ± 0.1) x 1CT 3 . 

The flavor bias of the detector Adet is due to asymmet- 
ric K + , K~ interaction cross-sections in the detector at 
low momenta. Such an asymmetry could produce a false 
CP asymmetry in the signal events. We perform a mea- 
surement of A^t in data using two independent meth- 
ods. The first approach determines this asymmetry from 
events in the toes sideband, 5.22 < toes < 5.27GeV/c 2 , 
which is dominated by continuum background with no 
expected CP asymmetry. The second approach uses a 
control sample where we replace the high-energy photon 
from the B decay with a high-energy 7r° with pcm > 
1.6 GeV/c. The same selection criteria used in the signal 
selection are applied, except for ir° and r\ veto require- 
ments, and the CP asymmetry in the toes sideband is 
measured. In both control samples we apply appropriate 
weights to the events to ensure that the fraction of each 
reconstructed final state is identical to that in the signal 
sample. We find the CP asymmetry measured using both 
of these approaches to be nearly identical, and average 
the two measurements to obtain Adet — — 0.007 ± 0.005. 
The mean value is used to shift the (N b - N^)/(N b + N-g) 



mean value, while the error contributes to the systemat- 
ics. The values of Adet computed in each X s mass region 
are reported in Table |TJ 

The shape of the BB and cross-feed background, deter- 
mined from MC, is also a potential source of flavor bias in 
the fit to the data. This background peaks broadly in the 
signal region, and a small shape difference as a function 
of flavor could create a false CP asymmetry in the sig- 
nal. We measure the size of this effect by correcting the 
BB and cross-feed shapes separately. The high-energy tt° 
control sample is used to study the uncertainty of the BB 
background shape. We use the differences found between 
the data and MC toes shapes in this control sample to 
correct the nominal BB background shape built from the 
MC. The biggest uncertainty in the cross-feed shape is 
due to the fact that JETSET does not reproduce the ob- 
served fragmentation structure of data. We thus correct 
the simulation shape using the fragmentation previously 
determined from BABAR data [2l[. We then construct 
new b and b binned PDFs using these corrected cross- 
feed and BB events and fit the data a second time with 
them. The difference between the nominal Aqp and Aqp 
from this fit, shown in Table []] is used as the systematic 
error from shape modeling of the B background. 

The systematic error arising from the continuum back- 
ground modeling is determined by varying the ARGUS 
shape parameters within the experimental errors, and is 
found to be 0.006 for the combined Mx s region. Sys- 
tematic errors due to possible differences in the signal 
shape between b and b events, CP content of the peaking 
background, and possible contaminations from b — > 
decays are all found to be negligible. Contributions from 
(D), AD and signal modeling are neglected due to their 
small impact on Ac p. The dominant systematic errors 
are therefore due to the uncertainties in the flavor bias 
of the detector and the background shapes as described 
above. 

The total systematic errors are calculated as the sum 
in quadrature of errors on Adet, systematic errors arising 
from the continuum, BB and cross-feed shape modeling. 
The results are shown in Table |TJ 

In summary, we measure the direct CP asymmetry in 
b -> s 7 to be A CP = -O.Off ±0.030±0.0f4 in the region 
0.6 < Mx s < 2.8 GeV/c 2 . This result represents the most 
accurate measurement of this quantity to date. The mea- 
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surement is consistent with zero CP asymmetry and with 
the SM prediction. The CP asymmetry in each Mx a re- 
gion considered in our study is also consistent with zero. 
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